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Abstract. It is widely believed that outflows of gamma-ray bursts are jetted. Some people also suggest that the 
jets may have structures like e(6) oc 9~ k . We test the possibility of X-ray flashes coming from such jets in this 
paper. Both qualitative and quantitative analyses have shown that this model can reproduce most observational 
features for both X-ray flashes and gamma-ray bursts. Using common parameters of gamma-ray bursts, we have 
done both uniform and nonuniform jets' numerical calculations for their fluxes, spectra and peak energies. It 
seems that nonuniform jets are more appropriated to these observational properties than uniform ones. We also 
give a spectrum and flux fit to the most significant X-ray flash, XRF971019 by our model. We also have shown 
that in our model the observational ratio of gamma-ray bursts to X-ray flashes is about several. 
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1. Introduction 

X-ray fiash(XRF) is a kind of recently identified explosion. 
Its most properties are qualitatively similar to those of 
gamma-ray burst (GRB) such as duration, temporal struc- 
ture, spectrum and spectral evolution except peak energy 
and flux. X-ray flash's peak energy and flux are lower, 
but their distributions just smoothly join the gamma-ray 
burst, there seems to be no obvious borderline between 
XRF and GRB. These similarities led to the suggestion 
that the X-ray flash is in fact "X-ray rich" gamma-ray 
burst (Kippen et al. 2003), maybe they have same origins 
except for different conditions. 

The similarity between XRF and GRB suggests that 
the X-ray flash might come from an off-axis nonuniform 
gamma-ray burst's jet(Woosley et al. 2003; Rossi et al. 
2002; Zhang & Meszaros 2002b). When a burst is observed 
at the center of the jet, it will be detected as a normal 
gamma-ray burst. But the burst tend to be "dirty" when 
it is observed at a large viewing angle (Zhang & Meszaros 
2002b), off-axis ejected matter takes less energy and has 
lower Lorentz factor. So its E p will shift to X-ray respon- 
sibility, and it will be observed as an X-ray flash. 

In this paper, we adopt a structured jet model where 
all the energy and mass of ejected matter per unit solid an- 
gle and the initial bulk Lorentz factor depend on the angle 
distance 9 from the center as power laws e(9) oc (9/9 c )~ k , 
m ej ((9) ex (6/0 c y k \ 7(6>) oc (9/9 c )~ kl (Meszaros et al. 
1998; Rossi et al. 2002). We take k = 2 for a nonuni- 
form jet, Rossi et al. have shown that 1.5 < k < 2.2 is the 



reasonable value for fitting observations well(Rossi et al. 
2002). 

Frail et al.(2001) had given the jet angles distribution 
with known redshift gamma-ray bursts. The jet's opening 
angle are range from 0.05 to 0.4 rad(Frail et al. 2001), 
and most common value is 0.12 rad(Perna et al. 2003). 
The gamma-ray energies released are narrowly clustered 
around £ 7 ~ 5 x 10 50 ergs(Frail et al. 2001). 

We introduce our model and give out some analytical 
solutions in Sect. 2. In Sect. 3 we present numerical results 
of spectra and fluxes for both uniform and nonuniform 
jets, and calculate the gamma-ray bursts to X-ray flashes 
observational ratio. Finally, we give a discussion and draw 
some conclusions in Sect. 4. 

2. The Model 

We consider a relativistic outflow where the energy per 
unit solid angle depend as power law on the angular dis- 
tance from the center 9 (Meszaros et al. 1998; Zhang & 
Meszaros 2002a; Rossi et al. 2002): 



e{9) 




(1) 



and the ejected matter per unit solid angle and the bulk 
Lorentz factor also depend on 9 as power laws: m c) {9) = 
m cj (0)(9/9 c )- k ^ 7 (0) = 7 (O)(#/0 c )- fcl (9 C < 8 < 9 S ). The 
deceleration radius at 9 is r d (9) = ( lm ,ff m ci ) 1/3 = 



r d (Q)(f )(-*+ 2fc O/3 . 
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All our calculations will be done at the time when an 
outflow just reaches its deceleration radius where the blast 
wave is formed. Because of the beaming effect of large 
Lorentz factor at this time, there is no obvious observa- 
tion difference between isotropic and anisotropic outflows. 
That means a jetted outflow with a viewing angle 9 V is 
observationally similar to an isotropic outflow with bulk 
Lorentz factor 7 = j(0 v ). So we can use the solutions from 
an isotropic explosion model (Sari & Piran 1999) to do an 
analysis by choosing different Lorentz factor at different 
viewing angle 9 V : 
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) 3 V (3)Fig. 1. The spectra of a non-uniform jet viewed from different 
viewing angles, here we adopt (9j = 0.1, 6 C = 0.02, k — 2, 
fei = k 2 = fc/2, r d = 4.0 x 10 16 cm, 7(0) = 500, p = 2.5, and 
the viewing angles 9 V change from to 0.16. 
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These equations describe the emission features from a 
shock between outflows and external mediums. Generally 
a external shock is not ideal for reproducing a highly 
variable burst(Sari et al,1998), but it can reproduce a 
burst with several peaks (Panaitescu & Meszaros 1998) 
and may therefore be applicable to the class of long, 
smooth bursts (Meszaros 1999). 

Electrons in the external mediums will be accelerated 
to a power law distribution. These electrons will product a 
broken power law spectrum with photon spectrum indexes 
ct(low) and /3(high) in the range —3/2 to —2/3 and — (p + 
l)/2 to — (p + 2)/2 through synchrotron emission (Katz 
1994; Cohen et al. 1997; Lloyd & Petrosian 2000). Here p 
is the power law index of accelerated electrons. There is 
no difference for an isotropic burst viewing from different 
direction, but for an anisotropic burst, E p changes from 
several KeV (or several eV, depends on 7(0), 6j, 9 C and 
k) to hundreds KeV (or several MeV) at different viewing 
angle, covers both GRBs and XRFs. 

From Eq.(2) and Eq.(4), we get: 
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Compared Eq.(5) and Eq.(6) we get 
F{6) ~ {vF v ) max {9) oc 



(5) 
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Here, 5 



_ fc+4fci 



only depends on the relation between k 



and k\. When k\ = k 2 = k/2 and 7(0) is a constant for 
every explosion, S — 3/2, that is the simplest solution. It 
will lead to the conclusion E p oc L 2 / 3 , which is close to the 
observational relation E p oc L 1 / 2 (Lloyd et al. 2000; Amati 
et al. 2002; Wei & Gao 2003). 
For an isotropic jet when 

(7(1 - 0008(0, - A<^-3 
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Fig. 2. The spectra of an uniform jet viewed from different 
viewing angles, here we adopt 8j — 0.1, rd = 4.0 x 10 16 cm, 
7 = 500, p = 2.5, and the viewing angles 8 V also change from 
to 0.16. Here lines represent to 0.08 are overlapped. 



6» j )))- 1 (Yamazaki et al. 2002; Granot et al. 2002), in this 
case we find S is about 4 and E p ex L 1 / 4 . 

Outflows with a lower Lorentz factor, which is called 
as a "dirty" fireball or a failed gamma-ray burst, may 
also produce a X-ray flash(Hcise et al. 2003, Huang et al. 
2002). It has the same spectrum and flux as we have just 
given out for a nonuniform jet. A dirty fireball will draw 
the same conclusion. It means that we cannot distinguish 
our model from a dirty fireball model just by a single X- 
ray flash. If simply assume that the bulk Lorentz factor 
T oc L 1 / 6 , we get E p oc T 4 oc L 2 / 3 . That means maybe 
we cannot distinguish nonuniform jet model from dirty 
fireball model even by statistical properties. 

Here we have to point out that Eq.(3) do not take cool- 
ing of electrons into account which may cause the F v ^ max 
to decrease 1 to 3 magnitude. We will take the cooling ef- 
fect into account in our numerical calculation in the next 
section. 
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3. Numerical result 

We have given out some useful conclusions using a simpli- 
fied analysis. But for more realistic, the observed flux at 
frequency v is the integral of equal arrival time surface of 
a jet: 



F v = 



Here P'(y') 
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f(x) is synchrotron radiation power 



at frequency v' from a single electron in the fireball co- 
moving frame (Rybicki & Lightman 1979). D = ("f(9) (1 — 
jScosO)) -1 is the Doppler factor translating from fireball 
co-moving frame to observer frame. is the angle between 
direction of outflow and line-of-sight. The electrons' dis- 
tributions can be written as(Dai et al. 1999): 
l.For 7 C < 7c, mim 
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Fig. 3. A spectrum and flux fit to XRF971019 using our cal- 

= 0.02, 70 = 500, 
4.2 x 10 26 cm ~ 



culation program. Here we chose 8j = 0.1, 



6 V = 0.08, 
136Mpc. 



5.4, r d = 4.0 x W Lb cm, d L 
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C\ = C 2 7 m f n + , 
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r{9) = 



The equal arrival time surface at 9 is: 

(3cT 



1 -/3cos9 
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Where N c = |r 3 ni is the total number of electrons per 
unit solid angle, equals to the number of protons in swept 
ISM. 7 C is the critical electron Lorentz factor above which 
synchrotron radiation is significant (Sari et al. 1998). 
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Our numerical results have been shown in Fig.l and 
Fig.2. We choose 9, = 0.1, 9 C = 0.02, 7o = 500, r d = 
4.0 x 10 16 cm, 4 = lx 10 28 cm, e c = 0.1, e B = 0.1, p = 2.5, 
k = 2 for a nonuniform jet and k = for an uniform one. 

It seems that for a nonuniform jet, the spectra and 
fluxes fit both GRBs and XRFs observations fairly well. 

- lMeV, 
these are 



Viewed from the center, E p is about lOOKeV 



and flux is about 10 
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3. For 7c > 7c 

N(j e ) = CaJc P 7e,min < 7c < 7c,r 



L iVc 



where 
C 5 = - 



-Nr. 



1-P 

-c,„i„ -7e,m P ax 

We assume the bulk Lorentz factor ^(9) keeps a con- 
stant before the outflows arrive at their deceleration radius 
r d . Here we choose k\ = k 2 which makes r d a constant at 
different 9, this assumption will make calculations very 
simple. But notice that for different 9, the time for out- 
flows arrive at rd is different. We calculate the flux when 
outflow reaches rd at viewing angle 9 V , that time is: 



T = 



(1 — /3cos(max[# v — 9j, 0]))r d 
0c 



1 j(9)~ 1 r(9y 1 n^ 1 (10) the typical values for GRBs(e.g. in Fig.l, the cases for 
6> v = 0,0.02,0.04). While when viewed from off-axis, E p 
is about lOKeV - lOOKeV, flux is about 10" 8 ~" 7 erg • 
cm _2 s _1 (e.g. in Fig.l, the cases for 9 V = 0.04,0.06), these 
are the typical values for XRFs. When 9 V = 0.08,0.1, E p 
is about a few KeV, the flux seems a little lower, but still 
can be detected if the source distance is not so large. 

But it seems that for an uniform jet, the flux from 
the jet edge, where XRFs are thought to be from in this 
model, are too low. In this case the spectra peak at 10- 
lOOKeV, the fluxes are about 10~ 13 — 10 crg-cm _2 s _1 (e.g. 
in Fig.2, the cases for 9 V = 0.105,0.11,0.12). It can only 
explain nearby XRFs, such as z < 0.2(Yamazaki et al. 
2002). But XRFs are more likely have cosmological ori- 
gins(Heise 2002). 

Fig. 3 is a fit to the most significant X-ray flash 
XRF971019 using our nonuniform jet model. Here we 
choose p = 5.4 because the high vF v power law index 
is about -1.7. This value seems a little unusual, but sim- 
ilar to the case of GRB980425 whose high-energy power- 
law photo index is 3.8 ± 0.7(Galama et al. 1998). XRF 
is defined as X-ray transients with duration less than 
1000s which is detected by WFC(on BeppoSAX) but not 
triggered GRBM(on BeppoSAX). Actually, this definition 
may lead to strong selection effects on those transients 
with high power law index. 



(12) 



(13) 
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We assume that every explosion has same luminosity 
and jet shape, luminosity at center is L„(0), and at view- 
ing angle 9 is L v (9) = L u (ti)(9 /9 c )~ 2kl . Depending on the 
viewing angle 9, a burst can be detected only at the dis- 
tance D{9): 



D(9) < D matx (6) = ( 



MQ)(e/gc)- 2fcl )1/2 



4ttF„ 



(15) 



Here F^mm is the threshold of a detector. So the numbers 
of GRBs(A^qrbs) and XRFs(7V X rf s ) are: 
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n is the number of bursts per unit volume, 9 CI is the crit- 
ical angle, if viewing angle is larger than it the explosion 
will be observed as XRFs. We assume the peak energy 
from the jet axis £? P;C = lMeV, here we divide XRFs 
from GRBs with E p — 90KcV, then the critical angle 



= ( 



-^p,c 

90KoV 



N, 



GRBs 



) 4k i9 c . We get: 



2.6 



XRFs 



(16) 



This result depends on the parameters we have cho- 
sen. Changing these parameters in a reasonable range will 
change this result a little but not much. On the other 
hand, these parameters can be restricted by observed ra- 
tio. 

4. Discussion and conclusion 

It is more likely a gamma-ray burst outflow is jetted. The 
energy distributions are more likely smoothly deducing 
with angle distance from the center than uniform in the 
cone. In this paper, we use a structured jet model which 
the energy per unit solid angle decreasing as e(9) oc 9~ k . 
We reproduce the key observational features of XRFs and 
GRBs. Our calculations are based on the external shock 
model which seems not ideal for reproducing GRBs with 
highly variable temporal structure, but still adapt to those 
bursts with smooth light curves. And our calculations 
should also adapt to the internal shock model. 

We have chosen 9 C — 0.02 and 7(0) = 500 in our nu- 
merical calculations. But the lower limit to 9 C is about 
l/7max ~ 10~ 3 , and the bulk Lorentz factor in the center 
have an maximum value to 7 max ~ 10 5 (Piran 1999; Rossi 
et al. 2002). Thus E p may be in excess of MeV when the 
axis just point towards us though the probability for this 
case is small. 

We take k = 2 for a nonuniform jet, Rossi et al. had 
shown that 1.5 < k < 2.2 is the reasonable value for fit- 
ting observations well(Rossi et al. 2002; Zhang & Meszaros 



2002a). For the cases of 1 < k < 3, the calculated spectra 
and fluxes, also the relation between E p and L are similar 
to that of k = 2. Because of beaming effect and shape 
of equal arrival time surface, the radiation mainly comes 
from outflows just pointing towards observer. We will ob- 
serve a similar spectrum and flux at a larger (smaller) view- 
ing angle when k is smaller(larger). But it will change our 
solutions of GRBs to XRFs observational ratio. We also 
assume ki = ki = 1 in the calculations through this pa- 
per, k\ 7^ &2 is a more complex case which may lead to 
very complicated calculation. 

We find that the observational ratio of GRBs to XRFs 
is about several in our model. It predicts that more XRFs 
or soft GRBs will be found in the future when more 
sensitive instruments are launched. Barraud et al. have 
presented 35 GRBs/XRFs spectra from HETE-2 whose 
band is 4-700KeV(Barraud et al. 2003), the numbers of 
bursts with E p higher and lower than 90KeV in these 35 
GRBs/XRFs are 24 and 11, similar to our calculated re- 
sult 2.6. 

The observation had shown a correlation E p oc 
L 1 ' 2 (Lloyd et al. 2000; Amati et al. 2002; Wei & Gao 
2003). This result is still not well explained except for some 
arguments based on some very simple assumption (Lloyd 
et al. 2000). We find in our model, for a single explo- 
sion viewed from different angles, E p oc L 2 ^ 3 , close to 
the observational relation. For a dirty fireball the relation 
between initial Lorentz factor and explosion energy is un- 
certain. If simply assuming that T oc i 1 / 6 , it leads to the 
same conclusion we have presented. 

Our model cannot be distinguished from a dirty fireball 
model just by a single X-ray flash, but for these two models 
the statistical properties such as the observational ratio of 
GRBs to XRFs are different. In addition, the structured 
jet model is also different from a dirty fireball in their 
afterglows. A very obvious feature for a structured jet is 
that when 9 V /9 C and k is large, there will be a prominent 
flattening in the afterglow light curve, and a very sharp 
break occurring at the time 7 ~ (9 V + 9 C )~ 1 after the 
flattening (Rossi et al. 2002; Wei & Jin 2003). We suggest 
that such features are more likely to be found in an X-ray 
flash afterglow rather than in a gamma-ray burst afterglow 
if an XRF does come from an off-axis nonuniform jet, 
because an XRF has a larger 9 V /9 C value than a GRB in 
this model. 

Orphan afterglows once caused great expectations to 
testify GRB collimation(Rhoads 1997). In a nonuniform 
jet, orphan afterglows may be generated in two ways. One 
way is when viewing angle is out of the jet edge(e.g. in 
Fig.l the cases V = 0.12,0.14,0.16, the fluxes will in- 
crease to detectable values at later time). In this case, 
the GRBs/XRFs are undetectable due to beaming effect 
but afterglows are detectable which are less beamed. The 
other probable way is that the Lorentz factor along the 
line of sight is sufficiently small that the peak energy E p 
is below the X-ray band. This case did not appear in our 
calculations, it will appear if we choose a smaller value of 
6> c /#j or a larger value of k\. 
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We have neglected the evolution of bulk Lorentz factor 
and lateral expansion of the jet which would make the cal- 
culations more complex and we think these effects are not 
very important before outflows arrive at their deceleration 
radius. Huang et al. have given out an overall evolution of 
jetted gamma-ray bursts in detail(Huang et al. 2000), we 
suggest that one should take all these effects into account 
for more realistic calculations. 
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